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Foreword 


The  physical  reasoning  which  led  to  the  conclusion  by  the  author  that 
particle  trajectories  in  a  rotating  field  are  stable  in  the  plane  of  rotation 
functions  as  an  introduction  to  the  reader  of  the  concept  of  plasm*1,  confinement 
in  a  rotating  field. 

The  magnetic  field's  relation  to  Maxwell's  equations  is  discussed  and 
the  induced  electric  field  is  presented,  whereupon  the  equations  of  motion  and 
their  exact  general  solutions  are  given. 

The  orbits  of  particles  are  derived  for  special  limiting  cases  on 
the  basis  of  the  general  solutions,  and  the  orbits  are  compared  with  the  analog 
computer  solutions  of  the  equations. 

It  is  shown  that  a  rotating  magnetic  field  alone  is  not  enough  to 
confine  a  plasmaj  but  that  the  rotating  field  used  in  conjuction  with  other 
types  of  stationary  fields  should  result  in  an  effective  plasma  containment 
device. 


Finally,  two  additional,  physically-possible  rotating  fields  are 
introduced  and  their  properties  are  discussed  briefly. 


I .  Orientation 


Interest  in  plasma  physics  has  grown  exponentially  during  the  last 
decade;  and  during  that  time,  many  confinement  devices  have  been  proposed  and 
some  tested,  but  none  were  sufficient  for  confining  that  all-important 
substance — the  thermonuclear  plasma.  All  of  the  confinement  devices  to  date 
have  had  a  common  fault:  particle  leakage . 

The  basic  method  for  confining  a  plasma,  which  is  considered  by  some 
to  be  the  only  possible  way,  is  to  trap  it  in  a  magnetic  field.  A  major  problem 
in  plasma  physics,  therefore,  is  the  design  of  a  magnetic  field  which  has 
minimal  particle  leakage. 

This  paper  describes  another  leaky  field  configuration,,  but  there 
exists  a  possible  application  of  the  proposed  device  such  that  the  leakage 
is  acceptable  (see  Section  VI). 


II.  Introduction  to  the  Rotating  Field  Concept. 

The  reasoning  which  led  to  the  conception  of  a  rotating  field  for 
confining  plasmas  proceeded  as  follows: 


* 


Confinement  is  achieved  in  the  plane  of  rotation  but  not  always  in  the 
direction  normal  to  the  plane  of  rotation. 
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(a) 


Figure  1  -  Stepwise  Rotating  Magnetic  Field. 


(b) 


A  particle  spirals  in  a  uniform  magnetic  field  as  in  Figure  1(a), 
Hence,  a  particle  would  escape  from  a  uniform  field  of  finite  dimensions. 
Suppose  now  the  field  in  Figure  1(a)  is  suddenly  replaced  by  the  perpendicular 
field  in  Figure  1(b).  The  original  component  of  °.xial  velocity  is  converted 
to  azimuthal  velocity  and  vice  versa,  so  the  particle  is  again  moving  along 
lines  of  force.  If  the  field-rotation  process  is  continued  in  this  stepwise 
manner,  some  of  the  particle  should  possess  orbits  as  shown  in  Figure  2(a). 


(a)  (b) 


Figure  2  -  Particle  Orbit  in  a  Stepwise  Rotating  Magnetic  Field. 

But  some  particles,  because  their  velocities  are  in  the  wrong 
direction  at  the  instant  the  field  changes,  will  possess  orbits  as  in 
Figure  2(b). 

This  clearly  reveals  the  fact  that  particles  tend  to  follow  lines 
of  force  as  they  rotate.  Consider,  now,  a  magnetic  field  rotating  with  a 
uniform  angular  frequency  u)Q .  If  o)Q  is  slightly  smaller  than  the  frequency 
of  rotation  co^  (Larmor  frequency)  of  the  particle  in  the  field,  then  the 
particle  obviously  spirals  as  in  Figure  3. 


-  Spiraling  Orbit  in  Rotating  Magnetic  Field  with  coQ  < 


Figure  3 
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But,  if  is  greater  than  a  particle  is  turned  by  the  field  quicker  than 
it  can  move  along  a  field  line  and  is  therefore  bounded  within  the  plane  of 
rotation. 


This,  therefore,  implies  the  following  confinement  conditions: 

(i)  coo  <  co^:  Particles  are  unconfined  in  the  plane  of  rotation 
(ii)  >cojj:  Particles  are  confined  in  the  plane  of  rotation. 

ci)q  =  co^  is  called,  hereafter,  the  critical  frequency. 

In  the  confined  state,  most  particles  are  however  unconfined  in  the 
direction  normal  to  the  plane  of  rotation.  This  may  be  realized  hy  considering 
a  particle  which  is  instantaneously  moving  upward  in  a  field  with  “o  >wt* 

In  one  xy -plane*  orbit ,  the  particle  will  not  have  time  to  complete  a  Larfcor 
orbit  (.cince  to  >ujj)i  and  hence,  after  one  period,  the  particle  is  still 
moving  upward.  Hence,  the  particle  is  unconfined  in  the  z-direction.  The 
exact  condition  for  confinement  in  the  z-direction  will  be  derived  later. 


III.  Fields 


The  rotating  magnetic  field,  illustrated  with  respect  to  an  xy- 
coordinate  system  in  Figure  4,  is  desired  to  have  the  form  of  Equation  (l)  for 


non-relativistic  rotation. 

B  =  B  [i  sin  co 
—  o  “  I 


+  i  cos  co  t]  (l) 

o 

The  magnitude  of  the  field  is  the 
constant  B  .  Hence,  Equation  (1) 
represents°a  uniform  rotating  field. 


FIGURE  4.  Coordinate  System. 


A  necessary  condition  for  B  to  be  a  physically  possible  field  is  for  it  to 
satisfy  the  wave  equation: 


,2b  i  is. 


C2  3tc 


t 

See  Figure  4, 


(2) 


i  mno  uaoiATarr 


-  4  - 


BBD-937 


Substitution  shows  that  Equation  (1)  is  not  a  solution  of  Equation  (2),  which 
means  that  a  uniform  rotating  field  cannot  exist.  However,  it  is  proven  in 
Appendix  1  that  Equation  (l)  represents  an  approximate  solution  of  Equation  (2), 
the  approximation  being  valid  for  small  frequency-distance  products. 

The  exact  expression  for  the  rotating  field  is  derived  in  Appendix  1 


B  =  B  [i  cos  ay  sin  co  t  +  j  cos  ax  cos  co  t] 
-  o  -  J  o  *  o 


where  a  =  —  . 


As  ay  and  ax  ■+  0,  Equation  (3)  reduces  to  Equation  (l). 

The  size  of  field  necessary  to  contain  a  one  hundred  million  degree 
plasma,  given  by  the  condition  that  the  magnetic  pressure 


P  =  B  /8n 
m 


exceeds  the  plasma  pressure 


P  =  nkT, 

P 

is  on  the  order  of  30kG.  It  is  presently  beyond  the  state  of  the  art  to 
construct  a  30kG  field  which  rotates  at  the  Larmor  frequency;  but,  as  discussed 
in  Section  VI,  large  rotating  fields  may  not  be  necessary.  It  is  proposed  that 
a  small  rotating  field  can  be  superimposed  on  a  large  stationary  field,  resulting 
in  a  practical  device. 

The  vector  potential, defined  by  the  equation 


is  given  by 


&  =  V  x  A  , 


A  =  B  [y  sin  to  t  -  x  cos  co  t]  k 
-o  '1  o  o  - 


for  the  approximate  field.  The  induced  electric  field  is  calculated  by 
substituting  Equation  (5)  in  J^raday's  law: 


E  =  -  l 

c  at 


(6) 
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The  resulting  expression  is: 


co  B 

E  - - 2-2  [y  cos  co  t  +  x 

—  c  o 

The  exact  expression  for  Inconsistent 
the  following: 


sin  co  t J  (7) 

o  • 

with  Equation  (3) ,  is  given  by 


E  =  -  Bq  [sin  o y  cos 


co  t  ♦  sin  ctx  sin  co  t)  k 
o  o  — 


(7a) 


Equation's  (l)  and  (7)  are  the  field  equations  which  enter  into  the  equations 
of  motion  derived  in  the  next  section.  The  simplified  field  expressions  lead 
to  differential  equations  which  have  exact  solutions,  and  they  are  fairly 
accurate  for  coq  -  co^. 

Both  field's  (1)  and  (3)  con  ''.in  another  approximation  which  is  not 
quite  correct.  It  was  assumed  that  v/c  «  1  at  every  point  in  space  where  v 
is  the  speed  of  the  field  with  respect  to  inertial  space,  which  is  obviously 
a  bad  assumption  for  high  frequencies  of  rotation  or  for  large  distances  from 
the  origin.  Calculation  shows,  however,  that  the  approximation  is  fair'  for 
co  /ov  <  10  and  Vx2-«-y2  <  10.  For  small  frequency-distance  products,  the  rela¬ 
tivistic  field  expressions  reduce  to  Equations  (1)  and  (7).  All  statements  in 
this  paragraph  are  proven  in  Appendix  2. 


IV.  Equations  of  Motion  of  a  Charged  Particle  in  ar  Electromagnetic  Field, 
(a)  Laboratory  Coordinates 

The  equation  of  motion  of  a  particle  with  charge  q  and  mass  m  in  a 
magnetic  field  B  and  an  electric  field  E  is  obtained  by  applying  the  Lorentz 
force  equation  to  Newton's  second  law.  The  vector  equation  of  motion  is: 

m|*=q  [E  +  ^VxB]  (8) 


Using  the  expression  for  E derived  in  the  previous  section,  the  component 
equations  are  immediately  written  as: 

mV  =  -  v  cos  co  t 
x  c  z  o 

m^y  =  sin  w  t 

j  c  z  o 

mV  =  -  22  (yo)  cos  co  t  +  xco  sin  u  t)  +  ^  (V  cos  co  t  -  V  sin  co  t) 

Z  CO  O  O  O  CX  OV  n 


(9) 
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where: 


x  =  V 


y  *  V 
*  « 


z  =  V 


(10) 


Equations  (9)  are  relatively  complicated  and  their  solutions  are  not  immediately 
apparent.  It  is  clear,  however,  that  the  sines  and  cosines  are  introduced  by 
the  rotating  field,  and  therefore  fransforming  to  a  rotating  reference  frame 
should  remove  these  factors  and  simplify  the  equations  considerably.  This 
indeed  happens  wider  such  a  transformation. 

(b)  Rotating  Coordinates 

A  rotating  reference  frame  is  chcosen,  as  illustrated  in  Figure  5, 
according  to  the  orthogonal  transformation: 


M 

sin  o>  t  -  cos  w  t 

= 

o  o 

\y/ 

i  cos  to  t  sin  a  t, 

'  o  O' 

\ 

where X  and  Y  are  rotating  coordinates. 


Y 


FIGURE  5. 


(11) 


The  transformed  equations  are: 


X=-aaY+w2X  ) 

o  o 

*  ■  ♦  a-o*  ♦  “o2*  *  “Lvz 


(12) 


/ 


Rotating 
X(0)  =  Y(0)  =  IV 
X(0)  =  Y(0)  =  IV 
K  =  20 

Time  Scale  1:1 


0  -  Independ.  Var , 


>  mw)  lAteworr 
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where  z  •  Z. 

(c)  Time  Scale  Transfortaation 

Equation  (12)  may  be  represented  more  simply,  parametrically,  if 
one  defines 


% 

{s- 

“L 


and  transforms  the  time  scale  according  to 


U5Qt  a  5  'l  , 


The  resulting  equations  in  rotating  coordinates  are 


dY  JS> 

2? 


,  „  dX  w2_  „  dZ 

25  dt  +  &  +  *  B 


J_  dY 
?  dr 


The  equations  in  laboratory  coordinates  with  transformed  time  scale  are: 


V  cos 
z 


V  sin 
z 


5  (y  cos  +  x  sin  C't)  +  V  cos  -  V  sin 

x  y  , 


(16) 


*mm>  mraa  lamcatocy 
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V.  Solution  of  the  Equations  of  Motion. 

The  equations  of  motion  will  be  solved  exactly,  and  solutions 
corresponding  to  particular  initial  conditions  will  be  given  as  calculated  by 
hand  or,  in  the  case  of  more  complicated  solutions,  as  calculated  from  the 
differential  equations  by  the  analog  computer. 

First,  solutions  will  be  obtained  for  the  simpler  equations  in 
rotating  coordinates.  Then,  the  solutions  will  be  transformed  back  to 
laboratory  coordinates. 

Equations  twelve  are  given  in  operator  notation  as  the  following 
matrix  equation: 

D2-*2 
-  2£D 
0 


The  solutions  are  of  the  form 
characteristic  determinant, 

D2-*2 
-  2£D 
0 

(a)  Critical  Frequency 

The  solutions  to  Equation  (l8)  provide  the  critical  frequency  for 
confinement  which  was  discussed  in  Section  II.  Equation  (18)  reduces  to  the 
sixth  order  polynomial  equation: 

D2  [D4  +  (2£2+l)  D2  +  S2}  =  0  (19) 


2£D 

2  2 
u  -  K 


I •y^ 

e  where  the  D's  are  solutions  to  the 


2CD 


-c2 


0 

-D 


D 


=  0 


(18) 


The  first  factor  yields  a  double  root  D  =  0  which  corresponds  to  solutions  of 
the  form  constant  and  t.  It  will  be  seen,  later,  that  an  expression  of  the 
fora  T  occurs  in  the  Z-oolution  only.  Hence,  the  stability  of  the  X  and  Y 
solutions  depends  on  the  second  factor  of  Equation  (19).  This  provides  the 
four  roots: 


D  = 


+  J 


/l  +  2S2  ±  V  l+8€‘ 


2 


(20) 
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Solutions  are  stable  when 


1  +  2^  +  Vl  +  8S2  >  0, 


which  implies  that  £  >  1.  This  is  equivalent  to  the  critical  frequency 
statement: 

*o  >WL 

(b)  Solutions  in  Rotating  Coordinates 


(i)  Stable  Solution 


Define 


'1  +  Zt  +  yi+8i 


i  +  2C  -  n+8i 


Then,  for  £  >  1,  solutions  are  linear  forms  which  include  only  the  expressions: 

+j“2T 

constant,  x,  e  and  e 

Substitution  of  the  linear  forms  into  the  equations  of  motion  determines 
the  relations  among  constants,  and,  by  using  Euler's  identity. 

+j<P 

e  -  cos  <p  +  j  sin  <p  ,  (24) 


the  general  solution  in  sine-cosine  form  is  obtained.  The  general  stable 
solution  in  rotating  coordinates  is  therefore 

2  fco.  2+S2-ll 


2  U.2+S2-l 

=  ^  “  B.  cos  (V+1a.)  ’ 

i=l  L  ^  i  - 

=  )  B.  sin  (ci>. x  +  a. )  +  C 

/  .  l  ii 


i=l 

2  b. 


=  aT  cos  (c°iT  +  ai^  “  ^CT  +  D 


Figure  6  shows  a  sample  plot  of  the  z-motion  and  the  corresponding  x  and  y 
motion  for  C  /  0.  Figure  7  shows  z-motion  and  the  corresponding  x  and  y 
motion  for  C  =  0. 


FIGURE  6 


«mn  mwa  umiatott 
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FIGURE  7 


i  -  iov 
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The  time  scale  in  the  figures  is  represented  by  the  variable  6  given 


e  = 


where  T  is  defined  by  Equation  (l4) . 

Equation  (29)  and  (31)  are  combined  to  give  an  explicit  expression 
for  drift  velocity  in  terms  of  initial  conditions: 


„  _  o  iqt  *o 

zD  =  ,  1 

X“F 


I  +  £ 

V  _  =  - r*°  in  real  time, 

ZD  .  1 


Calculation  on  the  basis  of  the  above  equation  shows  that  a  particle, 
after  a  collision,  will  escape  at  rates  on  the  order  of  lO^cm/sec. 

This  implies  that  a  rotating  magnetic  field  alone  is  unsuitable  for 
confining  a  thermonuclear  plasma;  but,  as  developed  later  in  Section  VI,  the 
rotating  field  will  have  important  applications  to  other  plasma  devices. 

Figures  5»  6,  and  7  clearly  illustrate  the  xy-stability. 


(ii)  Unstable  Solution 


Define  a.  =  jco, 

l  l 


(i  =  1,  2), 


then  the  a.'s  are  real  numbers,  and  the  rotating  solutions  are  linear  forms 
which  consist  in  linear  combinations  of  the  expression  constant,  T,  +a^T  and 
+a?T  e 

e“  .  The  linear  forms  can  be  expressed  in  terms  of  hyperbolic  functions,  time 

and  additive  constants. 

The  general  xy-unstable  solutions  are  given  by: 


X  =  Y  a.K.  sinh  (a.T  +  p  ) 
L  ii  i  i 

i=l 

Y  =  ^  cosh  (a/r  +  p^  + 


=  -V  —  sinh  (a.T  +  p. )  -  ^CT  +  D 
L  a.  i  ri 


(35) 
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(c)  Solution  in  Laboratory  Coordinates 

The  general  solution  in  laboratory  coordinates  is  derived  by 
applying  Equations  (11)  and  (14)  to  Equations  (25),  (26),  and  (27).  The 
result  is: 


k2+*2-i 

-  I  -Ssr  Bi 

i=l  L  1  J 


sin  cos  (co^+a^)  -  cos  sin  (co^r+a^)  -  c  cos  £t 


f(oc2+£2-l 


*  -  I  Bi  .in  sin  <V*ai>  *  -W- 

i=l  L  1  J 

z  =  ^  cos  (co/t+aj  -  +  D. 


cos  £x  cos  (cOjT+a^)  +  C  sin  £x 


As  £*♦«•,  the  coefficients  in  square  brackets  approach  1,  OQ,  and  (^-to  ) 
(£-«2)  l/VT  ,  which  gives  the  limiting  solution 


x  =  G1  sin  (^ij-t  +  y^) 
y  =  G2  cos  (^t  +  Y2)* 


Equations  (40)  and  (4l)  may  be  combined  to  give  the  equation  of  the  path  in 
the  xy-plane: 


2  2 


72  *  fa  -  cos  (Vr2>  *  s“2  'VV 


G1  °2  12 


The  above  is  the  equation  of  an  ellipse.  If  Y-Yo  =  “  and  G  =  G2,  the  orbit 
is  circular.  It  is  interesting  that  the  sense  or  the  orbital  angular  momentum 
vector  is  only  a  function  of  the  initial  condii  ons,  and  it  does  not  tend  to  line 
up  with  the  vector  of  the  rotating  field. 

The  path  in  inertial  space  corresponding  to  Figure  5  is  an  ellipse  with 

G^  =  56.6  cm  . 

G2  =  1  cm 

r1  =  -i#i’ 

y2  =  90° . 


Amro  mrrva  lAWJtArom 
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The  computer  solution  (Figure  9)  gave  corresponding  results  for  the 
first  half  orbit;  but  by  that  time,  the  ellipse  had  rotated  in  the  direction. 
This  rotation  was  removed  from  the  idealized  elliptical  orbit  when  it  was  assumed 
that  g+i  For  the  orbit  in  Figure  9t  K  was  20— considerably  less  than  infinity. 

The  effect  of  including  C  'Equations  (37)-(39))  in  the  solution  is  to 
superimpose  on  the  ellipse  a  small  amplitude  wiggle  with  frequency  5.  Since  C 
varies  as  £  ,  the  effect  is  negligible  at  high  frequency;  but  the  wiggle  was 
visible  on  computer  runs  for  £  <  5» 

The  revolution  frequencies  of  the  particle  were  checked  in  both 
coordinates  systems  on  the  analog  computer  and  were  found  in  excellent  agree¬ 
ment  with  those  predicted  in  the  general  solutions.  The  oscillation  frequency 
of  the  orbit  in  rotating  coordinates  also  checked.  These  important  frequencies 
in  radians  per  second  are  listed  here: 

Angular  frequency  of  revolution  in  lab  coordinates  =  “j/VF 

Angular  frequency  of  revolution  in  rotating  coordinates  »  co 

o 

Oscillation  frequency  of  orbit  in  rotating  coordinates  =  V2" 


VI.  Applications  of  the  Rotating  Field. 

(a)  Properties  of  the  Hotating  Field 

The  following  are  properties  of  the  rotating  magnetic  field,  revealed 
in  the  computer  solutions  of  the  equations  of  motion,  which  suggest  possible 
applications  of  rotating  magnetic  fields  to  other  existing  plasma  confinement 
devices: 


A.  All  particles  possess  stable  motion  in  planes  parallel  to  the  plane 
of  field  rotation  when  the  angular  frequency  of  field  rotation  o>  exceeds 
the  Larmor  frequency  of  the  particle 

B.  Linear  motion  of  a  particle  in  a  rotating  field  is  partially  trans¬ 
formed  into  azimuthal  motion. 

C.  Motion  in  the  direction  normal  to  the  plane  of  rotation  consists  in 
a  constant  drift  superimposed  on  a  stable  oscillatory  component. 

(b)  Suggested  Applications 

The  rotating  field  has  obvious  applications  to  the  following  three 
devices,  each  of  which  will  be  discussed  in  some  detail: 

A.  The  Cusped  Geometry 

B.  The  Magnetic  Mirror 

C.  The  Stellerator 
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(i)  The  Cusped  Geometry. 

Particle  leakage  in  the  cusped  geometry  occurs  at: 

1.  the  line  cusp. 

2.  the  points  cusps. 

There  seems  to  be  two  possible  schemes  for  applying  the  rotating 
field  to  the  cusped  geometry  in  order  to  prevent,  or  at  least  decrease, 
leakage : 

1.  The  superposition  of  a  rotating  field  locally  In  and  parallel  to  a 
region  which  includes  the  line  cusp,  and  the  superposition  of  a  rotating 
field  locally  near  each  point  cusp  in  planes  parallel  to  the  line  cusp. 

2.  The  superposition  of  a  rotating  field  parallel  to  the  line  cusp 
throughout  the  field  of  the  cusped  geometry.  The  point  cusps  world 
have  to  reflect  most  of  the  particles  that  the  rotating  field  forces 
normal  to  the  rotation  in  order  for  the  method  to  be  feasable.  This 
would  evidently  impose  a  restriction  on  the  ratio  of  the  rotating  field 
to  the  cusp  field. 

The  system  in  (1)  will  now  be  considered.  If  a  rotating  field  is 
applied  locally  in  a  region  which  includes  the  plane  of  the  line  cusp,  particles 
tending  to  escape  through  the  line  cusp  will  be  swept  azimuthally  against  the 
cusp  field  lines.  Both  the  reflective  properties  of  the  cusp  field  and  the 
confining  properties  of  the  rotating  field  will  confine  the  particles  within 
the  line  cusp.  Simultaneously,  the  drift  velocity  imparted  to  the  particles 
by  the  rotating  field  will  move  them  from  the  rotating  field  into  the  portion 
of  the  cusp  where  the  particles  are  in  a  confined  state  until  they  again 
approach  a  cusp. 

A  local  rotating  field,  it  is  here  proposed,  would  be  applied  locally 
also  at  a  value  of  Z**  between  the  point  cusp  and  the  plane  of  the  line  cusp, 
nearer  to  the  point  cusp.  Particles  which  are  escaping  along  lines  of  force 
toward  the  point  cusp  will  be  deflected  against  the  cusp  field  line  by  the 
rotating  field;  but,  because  of  the  3-drift,  will  sometimes  pass  beyond  the 
rotating  field  where  they  will  reflect  against  the  cusp  field  because  of  the 
circular  motion  which  they  acquired  in  the  rotating  field. 

Section  VII  (b)  discusses  a  method  of  producing  a  rotating  field 
which  is  especially  suitable  for  the  cusped  geometry  and  mirror  machines. 

(ii)  The  Magnetic  Mirror. 

The  rotating  field  could  be  applied  to  the  Magnetic  Mirror  for  the 
same  reason  application  to  the  point  cusp  of  the  cusped  geometry  seems  possible: 
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Particles  escaping  along  lines  of  force  are  given  azimuthal  motion 
by  the  rotating  field  which  causes  them  to  reflect  from  the  curved  field  lines, 

A  rotating  field  would  be  superimposed  upon  each  mirror  as  shown  in 

Figure  11. 


FIGURE  11.  The  Rotating  Magnetic  Field  Applied  Locally  to  the 
Magnetic  Mirror  Field. 


(iii)  The  Stellerator 

The  property  of  the  rotating  magnetic  field  which  makes  its  application 
to  the  stellerator  feasibile  is  the  following: 

Particles  have  stable  orbits  in  a  rotating  magnetic  field  when  u>q  >  a>^. 

The  rotating  magnetic  field  would  be  used  as  a  focusing  device  in  this 
application  to  return  particles  which  have  drifted  from  the  axis  of  the  stellerator 
because  of  the  curvature  and  field  gradient  drifts  which  are  inherent  in  the 
machine.  A  convenient  position  for  the  rotating  fields  would  be  on  opposite  legs 
of  the  stellerator  as  indicated  in  Figure  12. 


Stellerator  Field 


~>Rotating  Field 


FIGURE  12 


The  Rotating  Magnetic  Field  Applied  Locally  to  the  Stellerator 
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VII.  Alternate  Rotating  Fields. 

It  was  shown  in  Section  III  that  an  exactly  uniform  rotating  magnetic 
field  cannot  exist.  It  is  in  order,  therefore,  to  consider  physically  possible 
rotating  fields  and  to  note  a  few  of  their  properties. 

(a)  Crossed  Electromagnetic  Field 

A  rotating  magnetic  field  is  formed  by  (1)  the  superposition  of  a 
plane-polarized  EM  wave  which  progresses  parallel  to  the  Z-axis  upon  (2)  a 
plane-polarized  EM  wave  which  travels  parallel  to  the  y-axis  and  has  the 
right  phase  relation  with  the  previous  wave.  This  system  of  waves  will  be 
referred  to  as  a  "crossed  EM  field". 

The  correct  field  which  fulfills  the  arbitrality  imposed  coordinates 
defined  in  the  previous  paragraph  is  expressed  by  Equations  (43)  to  (46): 

Wave  (1) 


Wave  (2) 


Si 

=  B 

0 

sin 

(a>  t  -  oz)  i 
o 

(43) 

h 

=  E 

o 

sin 

(tot  -  oz)  i 
o  — 

(44) 

—2 

=  B 

0 

cos 

(tot  -  ay)  i 

o  — 

(45) 

=  E 

0 

cos 

(®t-  CTy)  k 

o  — 

(46) 

The  total  field  at  any  point  in  space  is  given  by 


B  = 


E  =  E^  + 


The  result  is  a  non-uniform  rotating  magnetic  field  in  the  xy-plane 
and  a  non-uniform  rotating  electric  field  in  the  yz-plane.  If,cr  is  small,  ^ 
the  field  is  approximately  uniform.  For  example,  if  =  3x10  ,  then  o  =  10"  , 

The  equations  of  motion  for  a  «  1  in  coordinates  which  rotate  with 
the  magnetic  field  were  derived  as: 

•  •  •  o  . 

X  +  2w  Y  -  (o  X  =  w,  c  sin  to  t  cos  to  t  (49) 

o  o  L  o  o 

Y  -  2»  X  -  «  2Y  +  to  v  =  c  sin2  to  t  (50) 

o  o  L  Z  L  o 


•  •  • 

Z-airY  +  aira>X  =  carc  cos  a>  t 
L  L  o  b  o 


(51) 
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The  solution  of  the  characteristic  determinant  (Equation  52)  of  the 
associated  homogeneous  equations  gives  the  result  that 


,2  2 

1  -  co 

2co  D 

0 

0 

-2co  D 

2 

D  -  co  ' 

o 

0 

CO  CO 

-cotD 

L  o 

L 

=  0 


(52) 


The  Z  solution  has  a  component  which  is  again  linear  in  time,  which  makes  one 
believe  that  Z-drift  must  be  a  fundamental  property  associated  with  rotating 
magnetic  fields.  The  determinant  also  provides  the  result  that  the  X  and  T 
and  solutions  are  stable  for  all  co  *,  unlike  the  previous  case,  and  the 
characteristic  frequencies  are 

=  1 +  ca^2  (53) 


co-  =  co  (54) 

d.  o 


The  inhomogeneous  equations  can  be  solved,  but  this  labor  has  not 
yet  been  attempted. 


The  characteristic  determinant  has 
there  are  solutions  of  the  form 


the  solution 


D 


=  0  which  implies  that 


x  =  A1  +  B^T 
y  =  A2  +  B2t 
2  =  Aj  +  Byt, 

but  substitution  shows  that  A^,  B^,  and  B2  are  all  zero. 
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(b)  Circularly-polarized  Electromagnetic  field. 


A  rotating  magnetic  field,  as  well  as  an  RE  field,  is  associated  with 
a  circularly-polarized  EH  field  which  traverses  the  Z  axis.  Such  a  field,  which 
is  self-consistent,  is 


B  =  -i  B  cos  (co  t  -  crs;)  +  j  B  sin  (cot  -  az) 

—  —  o  o  '**0 


(55) 


E  =  i  E  sin  (go  t  -  Oz)  +  j  E  cos  (cot  -<*z)  (56) 

—  —  o  o  “o 


The  equations  of  motion  in  laboratory  and  rotating  coordinates  are,  respectively, 

(57) 


V  B 

rax  =  qE  sin  (cot  -  az) - - —  sin  (cot  -  az) 

o  c 


7  B 

my  =  qE  cos  (cot  -  az) - cos  (cot  -  az) 

o  c 


V  B 
*•  x  ( 
mz  =  - —  sin 


V  B 


x  °  ~in  (cot  -  az)  +  °  cos  (cot  -  az) 

c  c 


and 


X  +  2ooY-co  X  =  -(e-  co.V„)  sin  aZ 
o  o  L  Z 


Y-2wX-co  Y=  (e  -  coV)  cos  aZ 
o  o  t  G 


Z  =  “  (Y  -  co  X)  -  (X  +  go  Y)  sin  aZ 

L  0  0 


(58) 

(59) 

(60) 

(61) 
(62) 


where  e  = 


qE 


m 


The  XY  rotation  in  the  Z  dimension, 


cos  oZ 
sin  aZ 
0 


sin  a  Z  0 
cos  a  Z  0 
0  1 


(65) 


removes  the  sine  and  cosine  coefficients  from  the  equations  to  give  Equations 
(64)  to  (66). 
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X  +  ouy  + 

U  -  o\y  - 
r  +  c^u  - 


sour 
2oir . 

°V* 


•  cr^Y^X  + 

•  + 

+  CO.  CO  X  : 


2co 

o 

2<o 

o 

.  0 


[-U+ctXy]  -  co  2X 
o 

[1+ouyJ  -  co  2u 

o 


0 

V-  • 


(64) 

(65) 

(66) 


No  obvious  result  is  immediate  except  that  for  small  a,  the  magnetic  field 
equations,  if  one  replaces  i  by  (-i),  reduce  to  (43)  and  (45).  Thus,  one 
expects  the  same  behavior  of  the  equations  of  motion  for  small  a  in  this  case 
as  in  the  crossed  field  case. 

If  interest  is  generated  in  regard  to  particle  confinement  in  crossed 
fields  or  circularly-polarized  fields,  it  may  be  feasible  to  attempt  to  solve 
the  equations  on  the  analog  computer  to  check  solutions  for  values  of  a  other 
than  the  limiting  one. 

It  should  be  noted  that  the  circularly-polarized  field  has  the  correct 
form  for  producing  a  rotating  field  in  planes  parallel  to  the  line  cusp  if  such 
a  wave  is  transmitted  along  the  axis  of  symmetry  of  a  cusped  geometry  machine. 


VIII.  Appendices . 

Appendix  1.  Derivation  of  the  Rotating  Field  which  Reduces  to  the  Given 
Uniform  Field  for  Small  a. 

The  equations: 

(Al) 
(A2) 


B  =  B  [i  cos  ay  sin  co  t  +  j  cos  ox  cos  co  t] 
—  o  —  o  o 


E  =  -B  [sin  ay  cos  co  t  +  sin  ox  sin  co  t]  k 
—  00  o  — 


are  derived  here  by  assuming  that: 


B  =  B  [i  sin  co  t  +  cos  co  t]  is  correct,  deriving  an  expres¬ 
sion  for  E  by  Faraday^s  law  and°obtaining  a°second  expression  for  B  by  Ampere's 
law.  The  total  field  is  then  equated  to  the  sums  of  the  first  and  second 
expressions  for  B,  and  the  iteration  process  is  repeated.  Additional  terras 
are  simultaneously  added  to  E. 


Faraday's  law 


and  Ampere's  law  are,  respectively, 


1  9B 

V  x  E  - - -  and 

—  -  c  art 


(A3) 


|  x  B 


1_ 

c  9ft 


(A4) 
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First,  assume  the  field  is  given  by 


B,  =  B  (i  sin  co  t  +  j  cos  to  t) 

— 1  O  ~  O'**  o 

and  substitute  B^  in  (A3),  Fhraday's  law,  to  obtain  the  electric  field 


E  =  -  -2-2  [y  cos  co  t  +  x  sin  (o  t]  • 
—  co  o 


Equation  (A6)  is  then  substituted  in  (A4),  which  gives  a  second  approximation 
B2  to  the  field  B: 

B  =  -  B  [iy2  sin  co  t  +  jx2  cos  co  t]  .  (A7) 

2c  0  0  ° 


Both  and  B.,  are  approximate  solutions  to  Maxwell's  equation  [(A3)  and  (A4)]  , 
as  B^  i  B.,  is  an  approximate  solution.  B  +  B.,  is  substituted  back  into  (A3) 
to  give  a  new  approximation  to  E: 

E_  =  -  B  (— y  -  )  sin  co  t  -  (— )  cos  co  t  k  (A8) 

■*  °  L\C  3.=3/  0  l'  3!C3/  »J- 

If  the  above  process  is  continued,  the  fields  are  represented  by 
the  infinite  series: 


r  /  22  44  6  6  \ 

•  -  v 


/  2  2  4  4  6  6  > 

+  i  +  ... 

V  2 ! c  4»C  6lc  y 


COS  co  t 
o 


/coy  33  55  77 

B  -ML  ,  Ml  _  “oL.  *  .. 


°  c  x'J  si-? 


3 ! c  51c'  7ic 


.  cos  co  t 
/  o 


3’c'  5!c5  7ic7/  ° 


sin  co  t  Ik 


(A10) 


The  expressions  in  parentheses  in  B  and  E  converge  to  sine  and  cosine 
respectively,  so  that  Equations  (Al)  and  (A2)  are  proven  to  be  correct. 
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Appendix  2.  Derivation  of  Relativistic  Fields. 

For  K  =  10,  the  velocity  of  the  magnetic  field  lines  approaches  the 
speed  of  light  at  about  100  cm  if  K  is  referenced  with  respect  to  conditions 
for  confinement  of  a  deuteron  at  100  million  degrees.  This  implies  that  the 
situation  must  be  described  by  relativistic  equations  for  accuracy. 

Maxwell’s  equations,  assumed  to  be  Lorentz  invariant  by  the  first 
postulate  of  special  relativity,  have  the  same  form  in  a  uniformly  translating 
frame  of  reference.  The  lorentz  field  transformations  are  given  by* 


B'l 

E’||  =  E|| 

E'  ,  =  — L=.  (E  +  i  Vx3) 

1  vt(F  "  c  1 


(All) 

(A12) 


(A14) 


Where  I  I  denotes  components  parallel  and  J.  components  perpendicular  to  the  axis  of 
translation,  and  p  =  V/c. 

With  the  assumption  that  most  of  the  motion  occurs  in  the  ^-direction 
equations  (All)  to  (Al4)  are  in  cylindrical  coordinates: 


B  ' 
r 


co  r 

-2-  (0xE) 

c  —  —  r 


Be'  =  Be 


B  ’ 
z 


co  r 

—  (9xE) 

c  —  —  z 


E  ’ 

r 


(0xB) 
- r 


Ee'  =  Ee 


e  • 

z 


(0xB)„ 


P 


co  r 


(A15) 

(Al6) 

(A17) 

(A18) 

(A19) 

(A20) 


where 
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For  simplicity,  the  low  velocity  fields  are  assumed  in  a  rotating  frame  of 
reference  since  it  is  known  the  low  velocity  limit  will  be  taken  anyway. 

In  order  to  have  a  uniform  rotating  field  in  inertial  space,  the 
field  in  rotating  coordinates  is  the  constant  B  and  the  electric  field  is 
zero.  0 


This  gives: 


E  '  =  0 
r 

V  «  o 

E  •  =  0 
z 


3  '  =  B  sin  (0  -  at) 
r  o  o 

Bq1  =  Bq  cos  (0  -  a^t) 


\ 


>  (A21) 

j 


/  (A 22) 


J 


Combining 


Equation  (A15)  to  (A22)  gives  the  results: 


B  cos  (0  -  at) 


(A23) 


B«  =  B  sin  (0  -  at)  (A24) 

*-  o  o 

a  r  B  cos  (0  -  a  t) 

E  =  _2 - 2 — - 2 —  (A25) 

z  cn-p* 


In  the  low  velocity  limit  limit  (A23)  to  (A25)  reduces  to  the  low 
velocity  limit  equations dn  Section  III. 

If  one  sets  £  =  10  and  r  =  10,  one  sees  that  P  =  .03  and 

which  implies  that  using  the  non-relativistic  equations  gives  fairly  good 
accuracy  up  to  this  limit  of  a  r. 
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Appendix  3*  Energy  Integral. 

If  the  first  equation  of  (12)  is  multiplied  by  X,  the  second 
multiplied  by  Y,  the  third  multiplied  by  2,  and  the  resulting  equations  are 
added  and  integrated;  the  following  conservation  law  results: 

2 

~  (X2  *  Y2  +  Z2)  -  ~  (X2  +  I2)  =  Constant. 


This  is  the  mathematical  statement  of  an  interesting  property  of  rotating 
fields:  the  kinetic  energy  with  respect  to  a  rotating  frame  of  reference 
minus  the  centrifugal  potential  is  constant. 
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